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Previoisly, we reported calculations of the high-temperature series for thermodynamical suscep- 
tibilities towards a number of possible ordered states in the t-J-V model. Due to an error in the 
calculation, the series for <i-wave supeconducting and extended s-wave superconducting orders were 
incorrect. We give the replacement figures. In agreement with our earlier findings, we still find no 
evidence of any strong enhancement of the superconducting susceptibility with decreasing temper- 
ature. However, because different Pade approximants diverge from each other at somewhat higher 
temperatures than we originally found, it is less clear what this implies concerning the presence or 
absence of high-temperature superconductivity in the t-J model. 



Previoisly we reported calculations of the high- 
temperature series for thermodynamical susceptibilities 
towards a number of possible ordered states in the t-J- 
V model. Due to an error in the calculation, the series 
for d-wave supeconducting (d-SC) and extended s-wave 
superconducting (s-SC) orders were incorrect. The first 
terms of the correct series are 



(i) _ (i) 
Xs-SC — Xd-SC 



(l + x)(l-3x) 
81n((l - x)/2x)' 



Here, (3 = 1/T is the inverse temperature, and the doping 
x is the average number of holes per site. The SC suscep- 
tibilities at small x — > are suppressed not quadratically 
but only logarithmically in x. Note also that the value 
x = 1/3 is not singular, as the zeros in the numerator 
and the denominator cancel. 

We performed the corrected calculations of the series 
for thermodynamical susceptibilities with an extra order 
in powers of j3 up to the term with (3 W . The coefficients 
of the corrected series are available upon request. Figs.Q] 
- [3] represent the analyzed results for the corrected series 
and should replace Figs. 2 and 3 in Ref. [J. 

For physically plausible parameters, t > J, the ten- 
dency of the superconducting susceptibility to decrease 
with decreasing temperature is not quite as prominent 
as in our original work. The increased variance between 
different Pades relative to our former results indicate 
that the series can only be trusted down to tempera- 
tures T ~ J, whereas for t — 2 J, the mean susceptibility 
reaches a maximum at T/J ranging from T/J ~ 1.2 at 
x = 6% to T/J ~ 1.3 at x = 26%. Fig. [3] shows that 
at t = 3J the maximum is shifted up to temperatures 
T ~ 2 J. We also note that the effect of the n.n. repul- 
sion V on the superconducting susceptibility is much less 
dramatic than we originally reported. 

As a result, although the corrected series still suggest 
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FIG. 1: (Color online) Superconducting (d-wave) susceptibil- 
ity for the "unphysical" case t = J/2 with V = (dashes) 
and V = J/4 (dot-dash) from the HTS to 1/T 10 (diagrams 
with up to Ne = 11 edges). Shading represents the standard 
deviation of nondefective Pade approximants as discussed in 
Ref. QJ. The corresponding curves for s-SC susceptibility (not 
shown) are similar but somewhat smaller at low temperatures. 
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FIG. 2: As in Fig. [T] but for t = 2 J. Only V = curves are 
shown; curves for V = J/4 (not shown) are similar but more 
suppressed at low temperatures. 
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FIG. 3: As in Fig. but for t = 3 J. 



the absence of high-temperature superconductivity in the 
t-J-V model at t > 2 J, this conclusion cannot be drawn 
as strongly as in our original Letter 

We should also note the papers by Koretsune and 
Ogata 0] and by Putikka and Luchini [^| which appeared 
after Ref. [H In these papers, the authors use a similar 
technique to study not the thermodynamical suscepti- 
bilities but the equal-time correlation functions for the 
t-J-V model with the special value V = J/4. They com- 



puted the high-temperature series for both the uniform 
and real-space equal-time pairing correlation functions 
up to 12th order. These authors derive an instantaneous 
superconducting correlation length which increases with 
decreasing temperature, but which never grows as large 
as one lattice constant in the accessible range of T/J. 
Nevertheless, they interpreted their results as an indica- 
tion that the t-J model has a superconducting ground 
state. In our opinion, the failure of the thermodynami- 
cal susceptibility we have computed to show any strong 
tendency to grow with decreasing temperature in the ac- 
cessible range, as well as the extremely short correlation 
lengths involved, call into question the strength of the 
conclusions drawn in Refs. 1, S 
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We analyze the high-temperature behavior of the susceptibilities towards a number of possible ordered states 
in the t-J-V model using the high-temperature series expansion. From all diagrams with up to ten edges, reliable 
results are obtained down to temperatures of order J, or (with some optimism) to J/2. In the unphysical regime, 
t < J, large superconducting susceptibilities are found, which moreover increase with decreasing temperatures, 
but for t > J, these susceptibilities are small and decreasing with decreasing temperature; this suggests that the 
t-J model does not support high-temperature superconductivity. We also find modest evidence of a tendency 
toward nematic and d-density wave orders. 



The discovery of high temperature (high-T c ) superconduc- 
tivity in the cuprate perovskites launched a renewed effort to 
develop an understanding of the physics of highly correlated 
electronic systems. It is clearly significant that superconduc- 
tivity arises in these materials upon doping a nearly ideal, spin 
1/2 antiferromagnetic insulating "parent" state. Indeed, there 
is a prominent school of thought|jl|] that holds that high-T c 
superconductivity is more or less inevitable in a doped quasi 
two dimensional (2D) antiferromagnet; for this reason enor- 
mous effort has been focused on studies of the t-J model 
[Eq. ([!])], as the simplest model of a doped antiferromag- 
net. However, while the existence of antiferromagnetic or- 
der in the spin 1/2 Heisenberg model (the zero doping limit 
of the t-J model) is well established^, it remains uncertain 
whether the 2D t-J model, by itself, supports high-T c super- 
conductivity. In addition to antiferromagnetism and super- 
conductivity, various other types of order[|3j ^, |(| ^] have 
been or may have been observed in the cuprates, including 
nematic [Q[ (or spontaneous breaking of the point-group sym- 
metry), charge-stripe, [Q spin-stripe[[ij and d-density[^| wave 
(also called staggered flux or orbital antiferromagnetic) order. 
It is thus interesting to determine which if any of these orders 
are generic features of a doped antiferromagnet. 

In this paper we report the results of an extensive high- 
temperature series (HTS) study of the susceptibilities of the 
2D t-J model toward various short-period orders. (Long- 
period stripe order cannot be readily studied using these meth- 
ods.) Naturally, the results obtained in this way are only reli- 
able at moderately high temperatures. However, correspond- 
ing to any low-temperature broken symmetry state there must 
be a susceptibility which diverges at the ordering transition; 
unless the transition is strongly first order, this susceptibility 
will be large, and an increasing function of decreasing temper- 
ature even at temperatures well above any ordering transition. 

Put another way, the HTS is sensitive to relatively short- 
distance physics (the range is determined by the order to 
which the series is computed). However, since the super- 
conducting coherence length in the cuprates is thought to be 
around two lattice constants, it is reasonable to expect that the 
10-12 terms we have computed in this series are sufficient to 



probe the physics of the model on length-scales relevant to su- 
perconductivity. It is important to stress that HTS[|8| is free of 
finite size effects that plague other computational techniques. 

Specifically, we have computed the HTS for antiferro- 
magnetic (AF), J-wave superconducting (d-SC), extended s- 
wave superconducting (s-SC), nematic (N), and orbital anti- 
ferromagnetic (d-DW) susceptibilities, as defined in Eq. (^). 
To extend the temperature range over which the results 
can be trusted, we have made use of standard methods of 
resummation,|^] which we describe explicitly below. In out- 
line, what we do is to construct a set of Pade approximants 
of related Euler-transformed series [see Eq. (|3|)] with different 
values of the parameter 0q, eliminate all "defective" mem- 
bers of this set, and then average over the remaining series. 
The variance in this average gives an estimate for associated 
errors; where the variance gets large we conclude that the re- 
sults can no longer be trusted. It is possible that by biasing the 
series, using additional information about the low temperature 
state obtained from other methods, one might be able to ex- 
tend the results to lower temperatures. However, without such 
additional information, it is our experience that using differ- 
ent prescriptions for resummation, or even adding a few ad- 
ditional terms to the series, does not significantly change the 
results or increase their range of validity. 

There have been previous high quality series studies [| 
pH] , |j"2| , pj| ] of this same model which inspired the present 
work, but they primarily focused on extrapolating the results 
to T = 0. What distinguishes the present study from these ear- 
lier studies is (a) we have obtained susceptibilities that were 
not previously computed[0 [l4j] and (b) we have contented 
ourselves with studying the high temperature behavior of the 
model. 

Our principal findings are as follows: 1. The results we ob- 
tain are reliable, without apology, for temperatures T > J, 
and are probably qualitatively correct down to T ~ J/2. 
However, we have not found any method of analyzing the se- 
ries that we trust to any lower temperatures. 2. In the "unphys- 
ical" range of parameters, t < J, with V = 0, the strongest 
incipient order is AF — the AF susceptibility, xaf, is large and 
shows a strong tendency to increase with decreasing temper- 
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FIG. 1: Temperature dependence of the AF susceptibility for t = 
J/2 (dashes, x = 1%, 16%, 31%) and t = 2J (dot-dashed lines, 
x = 1%, 6%) from the HTS to 1/T 11 (diagrams with N E < 10). 
Shading represents the standard deviation of non-defective Pade ap- 
proximants as discussed in text. 

ature, although this tendency gets gradually weaker with in- 
creasing doping x (Fig. pi). The SC susceptibilities are largest 
at values of doping x £ 16% (both d-SC shown in Fig. || 
and s-SC, not shown [|i5j; typically Xd-sc > Xs-sc)- How- 
ever, the inclusion of an additional nearest-neighbor (n.n.) 
repulsion, V = J/4, is sufficient to strongly suppress the 
pairing fluctuations. 3. In the physical range of parameters, 
t > J, where the bare interactions between electrons are truly 
repulsive, both Xd-sc (Fig. |) and Xs-sc (not shown[|l5[) 
are small and decreasing with decreasing temperature already 
for T < 2 J; the pairing fluctuations are further suppressed 
by the addition of small n.n. repulsion. Based on these ob- 
servations we conclude that the 2D t-J model with t > J 
probably does not support high temperature superconductiv- 
ity. 4. For t > J, commensurate AF fluctuations are mod- 
erate for x = 1% but are dramatically suppressed already at 
x > 6% (Fig. 0). At larger x, we find that the d-DW (Fig. @) 
and the nematic (Fig. ||) susceptibilities are both moderate, 
show a weak tendency to increase with decreasing tempera- 
ture, and remain virtually unaffected by the addition of a weak 
n.n. repulsion, V = J/4. 

We compute the high-temperature series for the 2D t-J 
model defined on the square lattice, 

H = ~y] Uj c\ a c ](J + J (S; • - -Uifi^ + V mrij , (1) 

where tij = t is the hopping matrix element, Ci a is the 
electron annihilation operator, Sj = (l/2)c ia T rTa iCi a i and 
rii = c\ a Ci a are on-site spin and charge operators, and the 
doubly-occupied sites are projected out. The canonical t-J 
model|Jl|] can be obtained from Eq. ([[]) by setting the n.n. re- 
pulsion V = 0, while the version of the model used in pre- 
vious high-temperature series studies [T(| [ill |l3|] can be 
obtained by setting V = J/4. All results presented in this 



X 




0.5 1 1.5 2 2.5 T 



FIG. 2: Superconducting (d-wave) susceptibility Xd-sc for t — J/2 
with V = (dashes) and V = J/4 (dot-dash) from the HTS to 1 /T 9 
(diagrams with Ne < 10). Shading as in Fig. |l|. Pairing fluctuations 
are strong at V — 0, but are strongly suppressed already at V = J/4. 
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FIG. 3: Same as Fig. pj for t = 2J. Pairing fluctuations are weak 
already at V — (dashes), and decrease further with introduction of 
weak n.n. repulsion (dot-dash). 
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FIG. 4: Staggered-flux susceptibility Xrf-DW for t = 2J, at V = 
(dashes) and V = J/4 (dot-dash) obtained from HTS to 1/T 9 
(N E < 10). Shading as in Fig. [IJ 
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FIG. 5: Temperature dependence of the nematic susceptibility xn for 
t = 2,1 and V = (dashes) and V = J/4 (dot-dash) obtained from 
HTS to 1/T 9 (iV B < 10). Shading as in Fig. |l[ The n.n. repulsion 
term V has virtually no effect on xn- 



paper refer to these two values of V. 

The susceptibilities (per unit site) \ can be expressed via 
the second derivatives of the free energy with respect to ap- 
propriately chosen perturbing parameters, or as the irreducible 
thermodynamic correlation functions of (properly projected) 
operators O, 



Xo 



= 1/T. (2) 



We take the staggered magnetization Oaf = Ylr elC * r Sr> 
Q = (tt, 7r), for AF ordering, the anisotropic part of 
the kinetic energy On = J2 r ( c l+icr c rcr 



"r+y a 



h.c.)/2 for nematic ordering, the staggered orbital currents 



DW 



C r+y<r Crcr 



h.C.) /2 for 



d-UW ordering [|l6), and the isotropic (anisotropic) part of 
the uniform pairing Osc = J2r (A r ,r+x ± A r , T+ $) /2 for 
s-SC (<i-SC) ordering, where the pair annihilation operator 
Aij = CtfCji + CjiCii- Respectively, the first coefficients 

of the corresponding series are Xaf ~ 00 ~ x )l^ Xn = 

Xd-DW = X O ~ x), and X ( Jl S c = x£-sc = P 3 ?- 

The definition of susceptibilities in terms of the derivatives 
of the free energy offers a convenient way for constructing the 
cluster expansion ||]. For each inequivalent lattice placement 
C of a given diagram (connected graph) with ue edges, the 
relevant traces (grouped by the number of particles) are com- 
puted using block-diagonal matrices of the cluster Hamilto- 
nian and the operator O. The traces are combined to produce 
the coefficients of the inverse-temperature expansion of the 
cluster susceptibility xo{0, V] C) and the thermodynamic po- 
tential 0,(0, y; C) dependent on the variable y = z/(l+2z) re- 
lated to fugacity z = e^. After the subcluster subtraction we 
obtain the irreducible weights of the cluster [whose expansion 
starts with O(0 nE ~ 1 ) or higher power of 0]. Then, combining 
the cluster weights for diagrams with up to Ne edges, we gen- 
erate the series exact in the thermodynamic limit to ~ Ne ~ 1 



((3 Ne+1 for AF order since the operator Oaf is defined on the 
vertices). As the last step, we perform the Legendre transfor- 
mation to obtain the series for the free energy F(x,0), and 
reexpress[|l0|](b) the obtained series for x m terms of the av- 
erage hole density x = 1 + (dCl/d/j^T- 

All obtained series for %, fi, and F were carefully com- 
pared with series computed analytically to 3 , and also with 
series to /3 5 obtained by a direct differentiation of the free 
energy expression generated with an independently-written 
Mathematica[|l7|| program. Among other consistency checks, 
we verified the cancellation of low-order terms in the irre- 
ducible weights. We have also compared the obtained free 
energy and related specific heat series with those for the t-J 
model with V = 1/4 from Refs. |9| [fjj and also the specific 
heat series at x = with the corresponding series for the 
Heisenberg model jl8[]. 

A standard way of extrapolating a power series in j3 is to 
construct a ratio of polynomials p n (0) / Qm(0) with match- 
ing expansion in powers of 0, referred to as an (n, m) Pade 
approximant. However, when only a few first terms of the 
series are known, and without a detailed knowledge of the 
structure of the singularities of the function, the procedure is 
plagued by spurious divergences. Indeed, the coefficients of 
the power series are only weakly modified if the numerator 
and the denominator of the fraction have close roots. The ac- 
curacy of the extrapolation would not be affected by cancel- 
lation of such factors, which amounts to using smaller n and 
m. Furthermore, it is not generally clear whether a series in 
would give better extrapolation than a series in a related vari- 
able = f(0). As a result, it has become standard practice 
to average^] over a large number of "non-defective" Pade ap- 
proximants generated for a family of functions f(0, 0q) with 
different 0q, and use the corresponding dispersion to estimate 
the errors. A specific difficulty of extrapolating the series for 
susceptibilities [compared to non-singular quantities such as 
n(k)] is that interesting susceptibilities can actually diverge 
at a critical temperature. Additionally, the peaks in x develop 
at relatively large spatial scale, meaning that they become pro- 
nounced only at sufficiently high orders. 

To generate the curves in Figs. 1-5, for each x(P) at a given 
set of parameters we constructed a large number of Pade ap- 
proximants (all n, m with n + m > (n + m) max — 4) for 
both the original series and a number of series in terms of the 
Euler- transformed variable [ 



(3) 



= 0/(00 + 0), 



using 0.3 < 0o < 10. We then eliminated as "defective" 
the approximants with positive real roots in denominators and 
numerators [susceptibilities are nonnegative, see Eq. (Q), and 
we do not expect an actual phase transition with a divergent x 
at such high temperatures], as well as the approximants with 
close numerator-denominator complex root pairs, and calcu- 
lated the averages and the corresponding rms deviations for 
different temperatures T = 1/0. 

As an independent method of analysis, we have also looked 
at the behavior of the untransformed and most closely bal- 
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anced (to close to n) Pade approximants with even to. When 
they are non-defective, these "best" approximants produce 
curves that look qualitatively similar to the average, and lie 
well within the shaded regions of the plots. We have also tried 
a number of variations on the averaging procedure, testing on 
the AF susceptibility xaf at x = 1%, where we expect a 
prominent peak at low temperatures. In particular, we found 
that the hyperbolic tangent transformation (3 = tanh(/3//?o) 
which was popular in previous HTS studies |Tl[ ], strongly 
suppresses all peaks, effectively flattening the extrapolated 
functions for all susceptibilities we computed. 

To summarize, our results indicate that, although in the 
"unphysical" region, t < J, the t-J model displays a sharp 
increase in the pairing fluctuations as the temperature goes 
down, this is no longer the case for t > J. Also, super- 
conducting fluctuations are strongly suppressed by introduc- 
tion of a small n.n. repulsion V. Thus we conclude that 
high-temperature superconductivity is probably not a generic 
feature ||20|] of doped AFs. Apart from xaf at very small dop- 
ing, none of the studied susceptibilities display remarkably 
strong fluctuations in the "physical" range t > J as the tem- 
perature goes down to T > J/2, but xn and Xd-DW do show 
a modest enhancement. 

In the future, we intend to continue studying the t-J- 
V model within the high-temperature series approach, for a 
wider range of orders in hope of identifying the most relevant 
fluctuations in the pseudogap region. We also intend to study 
the interaction of order parameters, by looking at various sus- 
ceptibilities in modified models where an ordering (e.g., hop- 
ping anisotropy) is imposed at the Hamiltonian level. Within 
this general approach, we also plan to study a range of related 
models, in particular an array of coupled t-J ladders. 
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